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PKS 1004+13: A HIGH INCLINATION, HIGHLY ABSORBED RADIO-LOUD QSO - THE FIRST 

RADIO-LOUD BAL QSO AT LOW REDSHIFT?i 

Beverley J. Wills^, W. N. Brandt^ & A. Laor'' 
Draft version February 1, 2008 

ABSTRACT 

The existence of BAL outflows in only radio-quiet QSOs was thought to be an important clue to 
mass ejection and the radio-loud - radio-quiet dichotomy. Recently a few radio-loud BAL QSOs have 
been discovered at high redshift. We present evidence that PKS 1004-1-13 is a radio-loud BAL QSO. 
It would be the first known at low-redshift (zom ~ 0.24), and one of the most radio luminous. For 
PKS 1004-1-13, there appear to be broad absorption troughs of OVI, NV, Si IV, and CIV, indicating 
high-ionization outflows up to ~10 000 km s~^. There are also two strong, broad [^ 550 km s^^), 
high ionization, associated absorption systems that show partial covering of the continuum source. The 
strong UV absorption we have detected suggests that the extreme soft- X-ray weakness of PKS 1004-1-13 
is primarily the result of absorption. The large radio-lobe dominance indicates BAL and associated gas 
at high inclinations to the central engine axis, perhaps in a line-of-sight that passes through an accretion 
disk wind. 

Subject headings: quasars: absorption lines — quasars: emission lines — quasars: individual 
(PKS 1004-^13) ~ UV: galaxies — X-rays: galaxies 



1. INTRODUCTION 

About 11% of radio-quiet QSOs show strong broad 
absorption troughs (BALs) displaced to the blue of the 
corresponding broad emission lines of OVIA1034, Lya, 
NVA1240, SiIVA1400, and CIVA1549, indicating high- 
ionization outflows up to 0.1c - 0.2c. Except for the ab- 
sorption, the ultraviolet spectra of BAL QSOs are very 
similar to non-BAL QSOs. This, together with limits on 
resonance photons scattered into the absorption troughs, 
and partial covering deduced from spectropolarimetry and 
absorption-line ratios, suggests that BAL outflows are 
common to all (radio-quiet) QSOs, but are recognized only 
when BAL material intercepts our line-of-sight. BALs ap- 
peared to be a unique property of radio-quiet QSOs. Seem- 
ingly this was an important clue to central engine physics: 
QSOs appear to produce either a highly coUimated, pow- 
erful, relativistic jet, or a sub-relativistic uncoUimated 
wind that accelerates BAL clouds, but not both (Stocke 
et al. 1992, Weymann 1997). Recently, with the advent 
of systematic identification of faint (mJy) radio sources at 
1.4GHz, several radio-loud BAL QSOs have been discov- 
ered, so the incidence of BALs no longer changes at the 
conventional limit separating radio-quiet and radio-loud 
QSOs (Brotherton et al. 1998a; Becker et al. 1999). Nev- 
ertheless, BAL properties probably still provide important 
clues to the production of powerful radio jets (Weymann 
1997), and the axis defined by powerful radio jets may be 
used to relate the geometry of BAL regions to the symme- 
try axis of the central engine. 

Another clue to the physics of mass ejection in BAL 



QSOs, and therefore perhaps in all QSOs, is the extreme 
weakness of their soft X-ray fluxes (0.2keV-2keV; Green 
et al. 1995; Green & Mathur 1996); even their more pen- 
etrating harder X-rays are extremely weak (2keV-10keV, 
Gallagher et al. 1999). In order to understand the nature 
of soft X-ray weak (SXW) QSOs, that is, those having 
Uox > 2 (3000A to 2keV), we have investigated the UV 
spectra of all PG QSOs with archival lUE and HST ul- 
traviolet spectra (Brandt et al. 1999). PKS 1004-^13 is a 
member of that sample, and is also the only radio-loud. 
X-ray-weak member. 

PKS 1004-1-13 (4C 13.41, PG1004-M30) is a V=15.2", 
QSO of low redshift, z = 0.240, thus Mb = -25. 6^ It is 
a classic powerful lobe-dominant radio source - 1.2 Jy at 
1.4GHz, or logLi.4GHz ^ 33.6erg s~^ Hz~^, with radio- 
loudness, logi?* — 2.32 (as defined by Sramek & Weed- 
man 1980). Elvis & Fabbiano (1984) noted its exception- 
ally weak X-ray emission, deriving a spectral index >2.01 
between 2500A and 2keV. They favored an explanation 
in terms of an exceptionally strong optical-ultraviolet Big 
Blue Bump, rather than absorption of the X-rays. They 
also suggested that a relatively weak EUV ionizing contin- 
uum could result in weak BLR emission. Here we present 
and discuss evidence that PKS 1004+130 is not only a 
powerful radio source, but also a BAL QSO. 

2. OBSERVATIONS & MEASUREMENTS 



We have retrieved and combined the NEWSIPS- 
calibrated lUE data from the short-wavelength camera, in- 
cluding our own lUE observation (SWP04942, SWP10850, 

^ Based partly on observations by the International Ultraviolet Explorer satellite, collected at the Villafranca Satellite Tracking Station of 
the European Space Agency. 
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SWP16824, SWP27517 with observation dates UT 1979 
April 15, 1980 Dec. 1, 1982 April 24, 1986 Jan. 12 re- 
spectively) , weighting these according to the square of the 
signal-to-noise ratio. All observations were at low disper- 
sion in the large aperture, giving a resolution of ~1100 km 
s~^ from 1170A to 1980A. Wavelength scales were checked 
using the geocoronal Lya line. We have applied an ap- 
proximate correction for lUE artifacts and the wings of 
the geocoronal Lya line, by subtracting the SWP artifact 
spectrum appropriate to point source spectral extraction, 
as presented by Crenshaw, Bruegman, & Norman (1990). 

By good fortune, the hne-of-sight to PKS 1004-^13 hap- 
pens to pass only 34' from the center of the Leo 1 dwarf 
spheroidal galaxy, so Bowen et al. (1997) used this QSO as 
a background probe of Leo I's halo. They obtained high- 
resolution (O.uApixel-i) HST GHRS spectra of PKS 
1004-hl3 from 1276A to 1562A. We have retrieved the 
HST archival spectra, applying the same wavelength shift 
as Bowen et al. to convert to a heliocentric velocity scale. 

Both the shape and flux density level of the GHRS spec- 
trum agree well with the lUE spectrum, suggesting the 
correctness of the flux-density calibration, despite the pos- 
sibility of a 1500 A (observed wavelength) bump artifact in 
the lUE SWP data (Kinney et al. 1991). 

3. RESULTS 

The HST and coadded lUE spectrum for PKS 1004-^13 
is shown as the lower spectrum in Figure 1. Beneath is 
a standard deviation spectrum derived from the scatter 
among the individual lUE spectra. The systemic red- 
shift Zcm =0.2401 is derived from our and Stockton & 
MacKenty's (1987) measured wavelength of the narrow 
[O HI] A5007 emission line. Based on this, above the spec- 
trum arc shown the expected wavelengths of redshifted 
broad emission lines. As noted by Kinney et al. (1991): 
"Based on this spectrum alone, 1004-1-130 would not be 
classified as a QSO because of the very weak Lya emission 
line." Even the expected strong C IV A1549 emission is not 
clearly present. To demonstrate the unusual weakness of 
the broad Lya emission, we note that its rest equivalent 
width EW = I7A is the smallest of 19 low-redshift, radio- 
loud, lobe-dominant QSOs from the HST archives. The 
mean for this sample is 140A. The next smallest EW(Lya) 
is 56A for 3C 288.1, for which strong associated absorption 
is clearly present (Wills et al. 1995). Below, we compare 
the UV spectrum of PKS 1004+13 with that of a typical 
QSO. 

Kinney et al. (1991) noted unresolved absorption fea- 
tures at I534A and 1919A that they attributed to N V and 
CIV doublets at a redshift of 0.24. The high-resolution, 
high signal-to-noise-ratio GHRS spectrum shows that the 
components of the N V AA1238,1242 doublet are further 
split and that there is corresponding strong absorption 
in OVIAA1032,1038 (Fig. 2a, b). These two absorption 
systems are also seen, but less clearly, in Lya (Fig. 2c). 
All these features are also marked beneath the spectrum 
of Fig. 1. These systems have heliocentric absorption 
redshifts of Za — 0.2364 and 0.2387, and are intrinsically 
broad - respectively ^600 km s^^ and 500 km s~^ total 

^For PKS 1004-1-13, Mb = -25.6, and core-dominance = -1.7 (§4). 
6A; [O III] /H/3 intensity ratio, 0.15; and H/3 asymmetry, 0.06 (Boroson 
FWHM(H/3), 6100 kms"!; EW ([OIII]), 17A; [O III] /H/3 intensity ratio 



width, representing outflows of 50±50 km s^^ (no out- 
flow) to 1200 km s^^. It is therefore not surprising that 
the corresponding high-ionization emission lines should be 
suppressed over this velocity range. There is no sign of 
low-ionization absorption at these redshifts in the broad 
]V[gIIA2798 emission line (Antonucci, Wills, unpublished). 
Both OVI doublets appear optically thick (t ° ^ 6), 

with absorbing gas covering only 70% to 80% of the con- 
tinuum source (using the method of Arav et al. 1999). 
Thus the absorption is intrinsic to the QSO. The result 
for the N V doublets are consistent, but with larger uncer- 
tainties. 

We also note the probable existence of broad absorp- 
tion troughs at even greater outflow velocities - up to 
~ 10 000 km s~^. The wavelength range corresponding 
to to —10 000 km s~^ is indicated by shading in Fig- 
ure 1, for OVIA1034, Lya, NVA1240, SiIVA1397, and 
CIVA1549. The CIV associated and broader absorption 
are visible on the lUE two-dimensional images for the 3 
long-exposure spectra. These BAL-like troughs suggest a 
reason why the Lya-N V A1240 broad emission is so weak: 
it has been suppressed largely by N V absorption. The Lya 
BAL would appear to be much weaker than for the higher 
ionization species. Above the PKS 1004-1-13 spectrum, we 
show for comparison a spectrum of the well-known BAL 
QSO, PC 0946+301 (Korista & Arav 1997), redshifted to 
align the broad absorption features as indicated by the 
shaded strips. Note the suppression of broad Lya emis- 
sion by the NV BAL for PG 0946+301. Higher quality 
UV spectroscopy would confirm the reality of BALs in 
PKS 1004+13. 

One way in principle to distinguish absorption from 
emission, and isolate the absorption spectrum of 
PKS 1004+13, is to divide by a 'normal' QSO spectrum. 
We and others (e.g.. Wills et al. 1999) have shown that 
the greatest spectrum-to-spectrum differences can be de- 
scribed by luminosity relationships (the Baldwin effect), 
dependence on Boroson and Green's (1992) optical Princi- 
pal Component 1 (PCI), and radio core-dominance (Baker 
& Hunstead 1995, Vestergaard 1998). Therefore we choose 
to divide by a QSO spectrum that is quite similar in lu- 
minosity, radio core-dominance, and in optical PCI prop- 
erties. We chose 3C263^ with EW(Lya) = 114A, and 
EW(CIV) = 86A, from a sample of radio-loud QSOs 
(Wills et al. 1995). Figure 3 shows the result, and serves to 
illustrate the approximate strength of the N V absorption. 
This technique would give an accurate absorption profile 
only if the template spectrum were an accurate match to 
the emission line spectrum of PKS 1004+13, and in the 
unlikely situation that the absorbing gas covered the con- 
tinuum and the whole of the emission line region. Note 
that if the broad absorption in PKS 1004+13 had been 
further blueshifted, as in many well-known BAL QSOs, it 
would have been more clearly seen against a well-defined 
continuum. 

4. DISCUSSION 

PKS 1004+13 shows the following characteristics typi- 
cal of QSOs with high-ionization BALs: (i) absorption to 

Optical PCI properties are: FWHM(H/3), 6300 km s"!; EW ([OIII]), 
& Green 1992). For 3C263, Mb = -27.0; radio core-dominance, -1.1; 
0.24; H/3 asymmetry —0.1. 
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high outflow velocities (^10 000 km s^^) in high- ionization 
lines of O VI, NV, CIV, and possibly Si IV, with a weak 
Lya broad absorption, and suppression of the broad Lya 
emission line, (ii) judging by the region of the N V 'trough' 
in the high-resolution GHRS spectrum, the troughs are 
smooth; probably the C IV seen only in low-resolution 
lUE spectra will not resolve out in higher resolution obser- 
vations, (iii) we measure a non-zero balnicity index, 850 
km s~^ (cf. Weymann et al. 1991), (iv) there are multi- 
ple (2), high-ionization associated absorption systems, as 
commonly go along with BAL absorption (e.g.. Ogle et al. 
1999), (v) it is extremely weak in soft X-rays, as in other 
BAL QSOs (Green et al. 1995; Green & Mathur 1996), 
(vi) there is scattered- light polarization (see below), and 
(vii) there is no sign of continuum reddening (Figure 7 of 
Elvis & Fabbiano 1984). 

Unlike any BAL QSOs reported in the literature^, 
PKS 1004+13 shows classical lobe-dominant radio struc- 
ture (Kellermann et al. 1994, Miley & Hartsuijker 1978), 
with core dominance, log (core flux/lobe flux) = —1.74 
(5GHz, rest frame). Thus, for the first time, it may be 
possible to relate BAL QSO properties to the axis of the 
central engine. To put the core-dominance in context we 
note that PKS 1004-1-13 is also a 4C radio source, selected 
at a low frequency, 178MIIz. We can compare the core- 
dominance with the distribution for 3CR sources, also 
selected at 178MHz, as this quantity is distributed sim- 
ilarly for all FRII radio sources in 3CR and 4C (Hoekstra, 
Barthel, & Hes 1997). Of 28 3CR QSOs with sufficient ra- 
dio imaging, only two have clearly smaller core-dominance 
- 3CR68.1, an extremely reddened and highly scattering- 
polarized QSO with strong associated absorption (Broth- 
erton et al. 1998b), and 3CR351 - with strong absorption 
by ionized gas in the X-rays and strong, high-ionization, 
associated absorption lines in the UV (Mathur et al. 
1994). Of 66 3CR radiogalaxies, half have smaller core- 
dominance. Thus, in unified schemes for radio sources, 
this indicates for PKS 1004-1-13 an inclination angle skim- 
ming the dusty torus, ~ 50°. Thus, the BAL region in 
PKS 1004+13 appears to lie at a large angle to the central 
engine's axis, and the BALs could arise in an equatorial 
accretion disk wind. This also ties in with pictures of the 
polarization scattering geometry in which BAL gas lies at 
low disk latitudes. 

About a third of BAL QSOs show polarization at the 
^ 2% level, but this is rare in UV-selected non-BAL QSOs. 
Polarization increases towards the blue as expected for 
scattered light (Schmidt & Hines 1999). PKS 1004+13 
shows continuum polarization at the ^ 2% level, increas- 
ing towards the blue (Antonucci et al. 1996; Stockman, 
Moore, & Angel 1984). The scattered-light geometry sug- 
gested by this polarization may be related to the partial 
covering of the continuum that we deduce for the OVI 
absorbing region; partial covering has been shown directly 
by changes in polarization across the absorption troughs 
of some radio-quiet BAL QSOs (Schmidt & Hines 1999). 



PKS 1004+13 could have provided an interesting test 
of scattering geometry for BAL QSOs, through a compar- 
ison of the polarization position angle with the position 
angle of the radio axis (117°). Unfortunately, the optical 
position angle is ill-defined at present. Antonucci et al. 
and Stockman et al. both find different and wavelength- 
dependent polarization position angles, suggesting vari- 
ability, and probably more than one polarization mech- 
anism. 

Is the discovery of BALs in PKS 1004+13 consistent 
with the 11% rate of occurrence among radio-quiet QSOs 
(Weymann 1997) or among radio-loud QSOs discovered in 
deep radio surveys? We can address this question by not- 
ing that Brandt et al. (1999) have searched for associated 
absorption lines in all available HST and lUE spectra with 
good signal-to-noise ratio, for all z < 0.5 QSOs in the PG 
optically-selected sample. PKS 1004+130 is the only one 
of 12 radio- loud QSOs to show probable BALs - consis- 
tent with the usual discovery rates. Are BAL properties a 
function of radio-loudness? If strong absorption near Zem 
is a common characteristic, this may explain why more 
radio-loud BAL QSOs have not been discovered. 

PKS 1004+13 is very bright in radio through UV wave- 
lengths, and it is nearby, enabling more detailed spectro- 
scopic follow-up and imaging at high spatial resolution. 
The low redshift also means that the broad hydrogen lines 
and narrow line region (NLR) emission lines are acces- 
sible in the optical and near infrared. High-quality UV 
spectroscopy is needed to establish whether PKS 1004+13 
is indeed a BAL QSO, and would, in any case, lead to 
important clues to physics and covering factor of the ab- 
sorbing outflows, the BLR, and the continuum that ionizes 
it. Imaging of the relatively bright radio nucleus, polar- 
ized (scattered) light and [O III] emission may lead to fur- 
ther clues to the geometry. PKS 1004+13 may therefore 
be a key AGN in understanding the relationship between 
radio-loudness and BALs, and the geometry of BAL QSO 
outflows. 
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'^ Recently, another FR II BAL QSO has been found in the FIRST survey (Becker, private communication). 
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Fig. 1. — Below: The observed ultraviolet spectrum of PKS 1004+13, derived from HST GHRS and lUE SWP data, shown on a rest- 
wavelength scale. The GHRS data have been smoothed (compare with Fig. 2). The standard deviation spectrum is shown, only for the lUE 
data. Above: The smoothed HST spectrum of the well-known BAL QSO, PC 0946+301 (Korista & Arav 1997), on a redshifted scale (top) 
that aligns its BALs with our suggested BALs for PKS 1004+13. The ordinates are in 10^^^ erg s^^cm^^ A^^ and 10~^^erg s~^ cm~^ 
A-l for PKS 1004+13 and PC 0946+301, respectively. The BAL velocity range, to -10,000 km s^^, is indicated by the shading. Note 
the overlapping troughs of Ly a and N V A1240. The expected wavelengths of emission-line peaks are indicated next to the labels (for 
PKS 1004+13, 2om = 0.2401; for PG 0946+301, Zem = 1.222). The expected positions of narrower absorption doublets are shown below the 
PKS 1004+13 data. Dots show the positions of Milky Way interstellar absorptions. The gap near 980A corresponds to geocoronal Lya. 
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Fig. 2. — The narrower absorption features in PKS 1004+13, illustrating the two redshift systems, and the widths of the absorption features 
for the wavelength regions of a) O VI, b) N V, and c) Lyo. The expected wavelength of the emission-line peak is indicated above the spectra. 
The long-wavelength edge of the absorption corresponds to zero rest-frame velocity for the longest wavelength pair of the doublets. For the 
ordinate scale, see Fig. 1. 
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Fig. 3. — Comparison of a smoothed, scaled spectrum of a matched radio-loud QSO 3C263 (see text) with the ultraviolet spectrum of 
PKS 1004+13. Both are corrected for Milky Way reddening. In the lower panel we show the division of the spectra, PKS 1004+13 by 3C 263. 



